home *** CD-ROM | disk | FTP | other *** search
/ Celestin Apprentice 2 / Apprentice-Release2.iso / Source Code / C / Applications / UIFlow 1.0.1 / UIFlow Source / VSet2.0 / temp / vtesthi.c < prev    next >
Encoding:
C/C++ Source or Header  |  1991-12-12  |  3.4 KB  |  102 lines  |  [TEXT/????]
  1. /*
  2.     SAMPLE C PROGRAM
  3.  
  4.     HDF VSET 2.1
  5.     Jason Ng NCSA MAY 1991
  6.  
  7.     This program demonstrates the use of the high-level VSET write routines.
  8.     It shows how data can be stored in vdatas, and how vgroups can be created
  9.     from vdatas, vgroups and other HDF elements.
  10.  
  11.     This example creates a file "vtesthi.hdf" that may be viewed using the
  12.     vset utility "vshow".
  13.  
  14.     DETAILS
  15.  
  16.     This example shows how pressure data  and color data can be stored as a 
  17.     vset.  To store pressure data, a field named "PRES", is defined of 
  18.     float type. The routine VHstoredata() stores the pressure values in a vdata.
  19.  
  20.     Color data comprises 3 components (red, green, blue). These can be stored as
  21.     3 different fields in 3 vdatas. But this example shows that they can be
  22.     treated as ONE field, and stored together as a compound field called "COLOR".
  23.     The number of components of a field is called its order.  In this case, 
  24.     "COLOR" has order=3. (whereas "PRES" above has order=1).
  25.  
  26.     The routine VHstoredatam() must be used to store values of a compound field.
  27.     This is similar to VHstoredata() but has an extra argument for the order
  28.     to be specified.
  29.  
  30.     Finally a vgroup is created, and the ids of the created vdatas are stored
  31.     in the vgroup. This effectively groups the vdatas together. This example
  32.     also shows that you can insert another vgroup, as well as a non-vset element
  33.     (in this case, some element with tag=7777 and ref=1) into a vgroup.
  34. */
  35.  
  36. #include "vg.h"
  37. #define printf 
  38.  
  39. #define fs "vtesthi.hdf"
  40. #define NP 100
  41. #define NC 60
  42. #define ORDERC 3 /* 3 color components: rgb */
  43.  
  44. main (ac,av)   int ac; char**av; {
  45.   int n = ac;
  46.  
  47.   printf("%s: tests high-level routines\n", av[0]);
  48.   printf("Creates 2 vdatas, a vgroup, and a non-vset element\n");
  49.   printf("then link them all into another vgroup\n");
  50.   printf("Vdata 1 contains an order-1 float field PRES\n");
  51.   printf("Vdata 2 contains an order-3  integer field COLOR\n\n");
  52.   printf("The non-vset element has tag=7777, ref=1\n");
  53.   doit();
  54.   printf("all done. created file %s\n", fs);
  55. }
  56.  
  57. /* ------------------------------------------------------------------ */
  58.  
  59. doit() {
  60. DF * f;
  61. float pvals[NP];
  62. int cvals[NC][ORDERC];
  63. int i,j;
  64.  
  65. int pid, cid; /* refs of vdatas */
  66. int eid; /* empty vgroup's ref  */
  67. int gid; /* vgroup's ref  */
  68. int tags[10], refs[10];
  69. char *CLS = "EXAMPLE";
  70.  
  71. /* --- generate data here --- */
  72.   for(i=0;i<NP;i++)  pvals[i] = 100.0  + i * 0.001;
  73.   for(i=0;i<NC;i++) for(j=0;j<3;j++)  cvals[i][j] =   i + j*100;
  74.  
  75.  
  76. /* ---- open a new file --- */
  77.  
  78. if (NULL==(f=DFopen(fs,DFACC_ALL,0))) { printf("open %s err\n",fs); exit(0); }
  79.  
  80. /* ---- create 2 vdatas --- */
  81. pid = VHstoredata (f, "PRES", (unsigned char*) pvals ,NP, LOCAL_FLOATTYPE, "pressure vals",CLS);
  82. if (pid  == -1) { printf(" VHstoredata store PRES err. "); }
  83.  
  84. cid  = VHstoredatam(f, "COLOR",(unsigned char*) cvals, NC, LOCAL_INTTYPE, "rgb colors", CLS, 3);
  85. if (cid  == -1) { printf(" VHstoredata store COLOR err. "); }
  86.  
  87. eid  = VHmakegroup (f, tags, refs, 0, "This is an EMPTY vgroup", CLS);
  88. if (eid  == -1) { printf(" VHmakegroup err\n"); }
  89.  
  90. /* --- create a new vgroup to store the 2 vdatas and the empty vgroup -- */
  91. tags[0] = VSDESCTAG;  refs[0] = pid;
  92. tags[1] = VSDESCTAG;  refs[1] = cid;
  93. tags[2] = VGDESCTAG;  refs[2] = eid;
  94. tags[3] = 7777 ;      refs[3] = 1;
  95.  
  96. gid  = VHmakegroup (f, tags, refs, 4, "here is a vset with 4 links", CLS);
  97. if (eid  == -1) { printf(" VHmakegroup err\n"); }
  98.  
  99. /* --- close the file --- */
  100. DFclose (f);
  101. }
  102.